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Development of chemically stable proton conductors for solid oxide fuel cells (SOFCs) will solve several 
issues, including cost associated with expensive inter- connectors, and long-term durability. Best known 
Y-doped BaCe0 3 (YBC) proton conductors-based SOFCs suffer from chemical stability under SOFC 
by-products including C0 2 and H 2 0. Here, for the first time, we report novel perovskite-type 
Ba 0 .5Sro.5Ceo.6Zro.2Gdo.iYo.i0 3 _5 by substituting Sr for Ba and co-substituting Gd + Zr for Ce in YBC that 
showed excellent chemical stability under SOFC by-products (e.g., C0 2 and H 2 0) and retained a high 
proton conductivity, key properties which were lacking since the discovery of YBCs. In situ and ex- situ 
powder X-ray diffraction and thermo-gravimetric analysis demonstrate superior structural stability of 
investigated perovskite under SOFC by-products. The electrical measurements reveal pure proton 
conductivity, as confirmed by an open circuit potential of 1.15 V for H 2 -air cell at 700 °C, and merits as 
electrolyte for H-SOFCs. 



The finite nature of oil resources and the ever-growing energy demand necessitate the finding of an alternative 
energy conversion technology that is highly efficient and free of greenhouse gas emissions 1 ' 2 . Solid oxide fuel 
cells (SOFCs) utilizing oxide-ion conductors (e.g., Y-doped Zr0 2 (YSZ)) due to higher efficiency (up to 
—80%), fuel flexibility, and combined heat power generation are being considered as alternative over conven- 
tional greenhouse gas emission system for stationary and mobile applications 3,4 . However, high operating tem- 
peratures, typically between 800-1000°C, results in material degradation, coking in case of direct hydrocarbon 
fuels and sulfur poisoning in Ni-based cermet electrodes. The operating temperature of SOFCs is commonly 
dictated by the choice of electrolytes; hence, efforts have been focused on intermediate temperature (IT) (400- 
700 °C) ceramic proton conductors to reap many benefits especially with economic metal interconnects 5-8 . 
Among the known electrolytes, aliovalent-doped BaCe0 3 (BCs) have demonstrated high proton conductivity 
(~ 10" 2 Son" 1 at 700°C), but, their poor chemical stability to SOFC by-products such as H 2 0 and C0 2 restricted 
them from being considered for proton conducting SOFCs 910 . 

Persistent efforts to improve the key features of BCs have shown that dopant with higher ionic size compared to 
Ce increases their proton conductivity, while higher electronegative dopants increase their chemical stability. 
Yttrium remains as one of the best candidates for doping Ce in BCs, whereas Yb and Pr co-doping exhibited 
mixed ionic and electronic conduction 1112 . A comparison of the ionic radii and electronegativity (see supple- 
mentary Table SI online) suggests that both Y and Gd seem to be ideal for Ce site doping in BCs, which is further 
supported by computational studies using a 'mean field approach' where Gd 3+ and Y 3+ have showed the lowest 
solution energy for doping in Ce site 13 . On the other hand, Sr-doping for Ba is proven to increase the phase 
stability under water vapor 14 . In contrast to BCs, BaZr0 3 -based proton conductors show appreciable chemical 
stability, but have poor sinter- ability and normally needs very high temperature sintering (> 1700°C) that makes 
them unsuitable for electrode supported SOFCs 15 " 17 . Here, we demonstrate a novel perovskite-type 
Bao.5Sr 05 Ceo.6Zr 02 Gdo.iY 0 !0 3 _s (here after called as perovskite I") with excellent chemical stability under 
C0 2 and water vapor at the desired IT-SOFC operating condition along with high proton conductivity that 
can be used in practical proton conducting SOFCs. 
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Figure 1 | Phase characterization and morphological studies, (a) PXRD 
pattern obtained for the as-prepared, Bao.5Sro.5Ceo.6Zro . 2 Gd 0 .iY 0 .i0 3 _5 
(perovskite I) revealing the formation of cubic perovskite structure with a 
very small impurity peak near 29 values of 30°, (b) TEM image of 
perovskite I (c) HR-TEM (d) SAED pattern of perovskite I with their 
corresponding crystal planes indexed again suggesting the formation of 
single phase and (e) SEM image of perovskite I showing grains in the range 
of 2-3 \i and revealing the dense nature of the pellet. 



Results 

Perovskite I of the nominal formula Ba 0 .5Sro.5Ceo.6Zro.2Gdo.iYo.i0 3 -§ 
was prepared by solid-state (ceramic) method at 1450°C for 24 h. 
Figure la shows the powder X-ray diffraction (PXRD) pattern of as- 
prepared sample. Figures lb- Id show the TEM, HR-TEM and 
selected area electron diffraction (SAED) patterns of the perovskite 
I, respectively. The crystallite size obtained from TEM image was 
found to be in the range of 200 nm. Further, crystallite sizes are 
narrowly distributed and are highly crystalline as seen from the ring 
intensities in SAED 18 . SEM image obtained on the polished pellet 
reveals the highly dense nature of perovskite I pellet and grains in 
2-5 microns range as shown in Figure le. 

Figures 2a and 2b show the results of PXRD and TGA measure- 
ments in various stability studies, respectively. FT-IR measurements 
on the perovskite I before and after C0 2 treatment at 800°C are 
shown in supplementary Fig. SI. PXRD pattern in Figure 2a after 
24 h water vapor treatment of perovskite I do not show any addi- 
tional peaks corresponding to the formation of Ba(OH) 2 . We 
extended the treatment time to 48 and 168 h, but no major change 
was observed in their PXRD pattern, suggesting significant chemical 
stability under water vapor (see supplementary Fig. S2 online). TGA 
curves obtained on these samples suggest a maximum weight loss of 
<1% up to 900°C after 24 h. However, exposure to water vapor for 
extended periods (for 48 and 168 h) show a steady increase in weight 
loss with exposure time (see supplementary Fig. S3 online). Despite 
the increased water uptake, PXRD do not reveal any additional peaks 
other than the perovskite-phase lines. FT-IR study on perovskite I 
before and after water vapor treatment studies show only a minor 
peak at 855 cm" 1 which could be due to Ba-O bond vibrational 
modes (see supplementary Fig. S4 online) 19 . 

The ratio between the PXRD diffraction intensities of (Ino/l2oo)> 
(I110/I211) and (I211/I200) remains constant after various stability 
measurements, further confirming the structural integrity of the per- 
ovskite I (see supplementary Table S2 online). An additional indica- 
tion of perovskite I chemical stability is its appearance before and 
after exposure to C0 2 and H 2 0 vapor as it retains its green color 
despite the harsh chemical stability tests (see supplementary Fig. S5 
online). To understand the role of C0 2 partial pressure on chemical 
stability, we passed a mixture of 1 : 1 volume mixture of C0 2 and N 2 
with 30% humidification for 140 h at 600°C. The resulting perovskite 
I's PXRD pattern does not reveal any additional diffraction peaks 
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Figure 2 | Comparison of chemical stability of Bao.sSro.sCeo.eZro^Gdn.iYo.iC^-g (perovskite I) under various environments, (a) PXRD pattern of 
perovskite I after exposure to C0 2 at 800°C and H 2 0 vapor at 90°C along with in situ high temperature PXRD patterns obtained under 10% C0 2 /N 2 and 
10% H 2 /N 2 at 800°C reveal its chemical stability and (b) TGA curve obtained for perovskite I under air, C0 2 and dry H 2 all resemble each other closely 
suggesting no significant change in the chemical composition. TGA of H 2 0 vapor treated sample shows drop in weight from temperatures as low as 100°C 
suggesting the incorporation of H 2 0 in the crystal lattice. 
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corresponding to BaC0 3 or Ba(OH) 2 formation (see supplementary 
Fig. S6 online). 

AC impedance plots obtained for the perovskite I under different 
environments at 600°C are given in Figure 3a, along with the fitting 
results. Equivalent circuits used for fitting is given in supplementary 
Fig. S7 and verification of the fitting results were done by Kramers 
and Kronig (KK) analysis (see supplementary Fig. S8 online) 20 ' 21 . 
Further, the calculated capacitance values from equivalent series fit- 
ting showed the presence of grain-boundary and bulk contributions 
at low temperatures, while at higher temperatures the grain-boundary 
effects disappear as observed from the removal of semi- circle corres- 
ponding to capacitance 10~ 8 F (see supplementary Table S3 online). 
Arrhenius plots for perovskite I and Gd-doped Ce0 2 (GDC) (one of 
the commercial IT-SOFC electrolytes) are shown in Figure 3b. The 
activation energy obtained under wet hydrogen atmosphere for per- 
ovskite I is 55 kjmol" 1 significantly lower than the 94 kjmol" 1 
obtained under air, confirming its proton conducting nature under 
wet hydrogen. Fuel cell polarization (hydrogen-air cell) plots 
obtained with the perovskite I as an electrolyte and Pt electrodes 
by passing humidified hydrogen and air are given in Figure 4a. The 
open circuit potential (OCP) observed at 700°C was found to be 
1.15 V which is close to Nernst potential for water formation. 

Discussion 

PXRD study (Figure la) suggests the formation of a single-phase 
perovskite-type structure in space group Pm3m (no.: 221) with a 
lattice constant (a) of 4.282(3) A. The FFT analysis of the HR- 
TEM image further specifies the formation of cubic crystal lattice 
(See supplementary Fig. S9 online). This significantly lower lattice 
constant in comparison with other rare earths and Zr co-doped BCs 
could be attributed to the incorporation of smaller Sr ( 1 1 8 pm) for Ba 
(135 pm). For example, Gd, Pr, Sm and Zr-codoped BCs showed 
lattice constants in the range of 4.387-4.330 A 22 ' 23 . The d-spacing 
obtained from the PXRD measurements are matching well with that 
of values obtained from SAED pattern (Figure Id) 18 . The density of 
the pellets was found to be about 95% as measured by Archimedes 
method. In comparison, Y-doped BaZr 0 85 Y 0 15 0 3 _§ prepared by 
solid state method was reported to have a density in the range of 
75-80% and hot pressing was required to obtain densities >95% 9 ' 24 
which make BZY unsuitable for scalable SOFCs 25 . 



Ex-situ PXRD pattern (Figure 2) obtained after exposure to pure 
C0 2 at 800 °C for 24 h reveals that the perovskite I retain the original 
perovskite-type structure with no additional peaks. Similarly, high 
temperature in-situ PXRD pattern obtained under 10% C0 2 /N 2 
doesn't show any second phase, confirming that the perovskite I 
remains robust against chemical reaction to C0 2 . TGA curve 
obtained under pure C0 2 resembles that of the TGA curve obtained 
in air, further supporting in-situ and ex-situ PXRD studies that C0 2 
do not react with perovskite I in the whole temperature range despite 
favorable thermodynamic conditions. However, FT-IR measure- 
ments showed a very minor peak due to carbonate. Free energy of 
BaC0 3 formation as a function of C0 2 for BaCe 1 _ x Zr x 0 3 suggest 
about 70 mol% of Zr doping is needed to eliminate the carbonate 
formation at 800°C 5 . Hence, it is fair to assume that doping of com- 
paratively higher electronegative elements in both A and B sites of 
BCs have increased their resistance toward the acidic C0 2 and 
reduced the kinetics of BaC0 3 formation significantly. While proton 
conducting SOFCs may operate normally on external reformed 
hydrogen that produces lower C0 2 partial pressures, we carried 
out these stability measurements under extreme conditions to prove 
the superior stability of the perovskite I (Figure 2). Similar studies on 
YBCs reported to form carbonates at 500°C under C0 2 exposure as 
low as 2% 26 . It is important to mention that the flow rate of gases and 
the final sintering temperature of the investigated perovskite seem to 
be crucial for the observed chemical stability under C0 2 . 

The Gibbs free energy for the reaction between BCs and H 2 0 
vapor reaches the value of zero at 403°C and accordingly aliovalent 
doped BCs are reported to be stable at elevated temperatures 
(>500°C) under water vapour, while stability at low temperatures 
remained un- achievable 27 ' 28 . Our recent work on Sm and Gd + Pr co- 
doped with Zr in BCs could not achieve chemical stability under 
water vapor at 90°C even after 30% Zr-substitution for Ce 22,23 . But 
in the case of perovskite I, we could achieve water vapor stability even 
at low temperature of 90°C for extended duration (Figures 2a and b, 
see supplementary Figs. S2 and S3 online). For 168 h water vapor 
exposed perovskite I sample, the TGA showed increase in weight loss 
when compared to 24 h and 48 h water vapor exposed samples (see 
supplementary Fig. S3 online). 

Acceptor- doping in BCs can generate oxide ion vacancy according 
to the substitution mechanism: 
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Figure 3 | Conductivity measurement of Ba 0 .5Sro.5Ceo.6Zro.2Gdo.iYo.i0 3 _5 (perovskite I), (a) Nyquist plot obtained for perovskite I under different 
environments revealing the increase in resistance in the order of H 2 + 3% H 2 0 < air + 3% H 2 0 < air < N 2 + 3% H 2 0 < N 2 + 3% D 2 0 at 600°C. The 
lines passing through the data points are fitting results, (b) Arrhenius plot of perovskite I in various environments along with that of GDC measured in air 
revealing increased conductivity for perovskite I against GDC. Activation energy of 55 kjmol -1 is observed under H 2 + 3% H 2 0 environments. 
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Figure 4 | Polarization studies on Bao.5Sr 0 .5Ceo.6Zro.2Gdo.iY 01 03_5 (Perovskite I). Fuel cell polarization plot obtained using Pt paste as the cathode 
and anode and Perovskite I as electrolyte under humidified H 2 and dry 0 2 as the fuel and oxidant at 700°C and (b) Impedance Nyquist plot 
obtained for this cell under above mentioned conditions in the frequency range of 1 MHz to 0.1 Hz. The high resistance observed in the Nyquist plot is 
due to thicker electrolyte (~3 mm). 



Y-doping for Ce 

B4a(s) + 3 °o(s) > B<(s) + (1 -*)C&M 

+ *YU) + ( 3 " © ) °5« + © V "o(s) + © °*) 

where every 1 mole of M 3+ doping creates 0.5 mole of oxygen vacan- 
cies. On complete filling of these vacancies by H 2 0, TGA will result 
in a weight loss of —0.6%. However, 168 h exposure to H 2 0 vapor at 
90°C has resulted in weight loss of 3.5% which suggests excess water 
maybe adsorbed at surfaces or possibly at the grain-boundaries. SEM 
images of powdered samples after 168 h exposure are shown in Fig. 
S10 (see supporting information online). A similar study on 
BaCeo.g4Zro.01Smo.15O3- 5 resulted in a weight loss of >8% and in 
FT-IR, formation of peaks corresponding to Ba(OH) 2 in 24 h (see 
supplementary Fig. Sll online) 22 . Sm-doped BCs prepared in our 
group showed complete decomposition of the powders into amorph- 
ous materials in PXRD 22 . Thus, the increased water uptake upon 
extended exposure must be due to incorporation of water molecules 
in the vacancies and adsorption rather than formation of Ba(OH) 2 as 
also confirmed by PXRD measurements. The diffusivity of water in 
doped perovskites is reported to be slow at low temperature, which 
could be the reason for the slight increase in weight due to water 
uptake for 168 h compared to 24 h and 48 h (see supplementary Fig. 
S3 online) 2930 . 

Chemical stability under H 2 at elevated temperature is another 
significant factor as Ce based electrolyte materials tend to react under 
this highly reducing condition and result in lower cell potentials. HT- 
PXRD under 10% H 2 /N 2 further proves perovskite Fs possible utility 
as electrolyte as original perovskite structure is retained at 800 °C. 
Similarly, TGA curves under pure H 2 do not reveal significant weight 
loss and further proves the structural stability of the perovskite I 
under reducing conditions. Nevertheless, 0.13% weight loss in excess 
is observed for TGA in H 2 in comparison to that obtained in air at 
800°C which would correspond to —4.5% conversion of Ce 4+ to 
Ce 3+ . However, at the operating temperature of proton conducting 
SOFCs typically at 600°C, this conversion will be less than 2%. (see 
supplementary information for the calculation). Thus, based on the 
above results, we could conclude that perovskite I is considered to be 
highly stable under SOFC operating conditions, which is critical for 
its successful operation. 

Perovskite 1 shows the highest conductivity under H 2 + 3% H 2 0 
(Figs. 3a and 3b), while lowest under N 2 + 3% D 2 0 among investi- 
gated environments. In air, the conductivity is significantly higher 
under wet air than dry air suggesting proton conduction under wet 



conditions. Proton conductivity in doped BCs arises from the filling 
up of hydroxyl ions in the oxide vacancies (V 0 ): 

VZ {s )+H20 ig) +(% is) ~20H- 0{s) (2) 

At low temperatures, the conductivity is dominated by proton trans- 
port due to higher H 2 0 uptake. However, as indicated by TGA curves 
at high temperatures the protons leave the lattice thereby reducing 
the concentration of hydroxyl protons and contribution from oxide 
ion increases. Comparison of perovskite Fs conductivity against that 
of oxide ion conducting Gd-doped Ce0 2 (GDC) made from com- 
mercial powders (GDC 10 TC - Fuel Cell Materials) show two orders 
of magnitude increased conductivity at low temperatures. Never- 
theless, GDC is reported to show higher conductivity value in the 
literature 11 . A maximum ionic conductivity of 0.4 X 10" 2 Son" 1 is 
observed at 600 °C for perovskite I, where the in-house prepared 
GDC showed a conductivity of 0.3 X 10" 3 Son" 1 under similar 
experimental condition. BaZr 0 8 Y 0 2 O 3 _§ showed a conductivity of 
7.9 X 10~ 3 Son -1 at this temperature after hot pressing to ensure 
densification 24 . A comparison between the reported conductivity 
values of state of the art YSZ and perovskite I also show the superior 
conductivity at desired fuel cell operating temperature of below 
700°C U . 

At higher temperatures, the difference in conductivity obtained for 
perovskite I between wet H 2 and other atmosphere is diminished 
and could be attributed to proton dissolution from the crystal matrix 
at elevated temperatures 29 . An increased activation energy of 
77 kjmol -1 observed in 3% humidified air against the 55 kjmol" 1 
under wet hydrogen suggests possible mixed (protons and oxide 
ions) ionic conductivity. Similar study on acceptor-doped ortho- 
niobates showed that at temperatures higher than 700 °C, contri- 
bution of protons towards total conductivity decreased as observed 
by the decrease in proton transport numbers 31 . 

The open circuit potential (OCP) observed at 700 °C was found to 
be 1.15 V, which clearly reveals the highly dense nature of the mem- 
brane. Furthermore, it also confirms that the perovskite I is a pure 
ionic conductor. Proton transport number studies on similar per- 
ovskites demonstrated pure proton conduction up to 700 °C under 
wet conditions 1131 and we believe the same for the investigated per- 
ovskite I. Further, the small amount of Ce 4+ /Ce 3+ conversion 
observed under TGA has not induced significant electronic conduc- 
tivity as otherwise the OCP would have reduced due to electronic 
short circuiting. A maximum power density of 18 mWcm" 2 is 
achieved at 700°C for un-optimized thick perovskite I electrolyte 
with the use of Pt paste as electrodes. The Nyquist plots obtained 
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with perovskite I in fuel cell mode by passing wet H 2 and air are in 
accordance with the plots obtained in wet H 2 (Figure 3). A maximum 
conductivity of about 1 X 10~ 2 Son" 1 is achieved at 700°C under fuel 
cell operating conditions. 

Methods 

Preparation of Bao.5Sro.5Ceo.6Zro . 2 Gd 0 .iY 0 i0 3 _ s (Perovskite I). High pure 
(>99.9%) oxide and carbonate precursors, barium carbonate, strontium carbonate, 
cerium oxide, zirconium oxide, gadolinium oxide and yttrium oxide (Sigma Aldrich) 
were mixed in the appropriate ratio and ball milled for 6 h using 2-propanol as the 
solvent. The mixture was than dried and calcined at 1050°C for 24 h followed by 
further ball milling for 6 h. Thus, prepared powders were uni-axially pressed into 
pellets using an isotactic press at 200 MPa for 3 minutes and sintered with their 
parent powders for 24 h at 1450°C. 

Chemical stability measurements. Stability under water vapor at 90°C was 
performed by hanging perovskite I powder above water in a round bottom flask fitted 
with a condenser for various time durations. Stability under C0 2 gas was measured by 
placing Perovskite I powder inside a quartz tube and flowing pure C0 2 at a flow rate of 
10 seem for 24 h at 800°C. To simulate the SOFC operating conditions, 1 : 1 ratio of 
N 2 and C0 2 was purged through water that was kept at 80°C. PXRD and 
thermogravimetric (TGA) analysis (Mettler Toledo, TGA/DSC/HT 1 600) were 
performed on these samples before and after the stability tests. Heating rate during 
TGA measurements was 10°C/min. For in-situ TGA under C0 2 environment, the 
heating rate was decreased to 5°C/min and a flow of pure C0 2 was maintained. In-situ 
high temperature (HT) - PXRD measurements under various gas environments were 
performed in a high temperature reactor chamber (Anton Paar XRK 900) from 20 10° 
and 80° at a count rate of 3 s per step of 0.05°. FT-IR measurements were carried out 
on a Varian 7000 FT-IR spectrometer. 

Conductivity measurements. Sintered pellets were cut into smaller disks and both 
sides were grounded and polished to obtain the desired thickness and flat and parallel 
surface to one another. Pt paste was brushed on both side as electrodes and fired at 
800°C for 2 h. Conductivity was measured using AC impedance analyzer (Solartron 
electrochemical impedance spectroscopy; SI 1260) at various temperatures and under 
various environments. Prior to measurement, samples were held at the temperature 
of measurement for a minimum of 2 h and maximum of overnight. Single fuel cell 
polarization measurements were carried out by passing humidified H 2 and air on the 
anode and cathode sides, respectively. During heating 5% H 2 was supplied and is 
increased to pure H 2 after reaching the desired temperature. 
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